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ABSTRACT 

There are only about 65 R Coronae Borealis stars known in our Galaxy, and none in globular clusters. As these stars are thought to 
result from the merger of two white dwarfs, one would expect the higher stellar density of globular clusters to favor their formation. 
We have searched for such stars in Galactic globular clusters, as their presence in a specific category of clusters might provide more 
clues as to their formation. 

We selected from the WISE all-Sky source catalog all the stars within the tidal radius of the 150 globular clusters within 50 kpc, 
which is the distance to which RCB stars are detectable by WISE. The total number of stars selected in this way was 635989. We then 
successively applied the eight selection criteria of Tisserand (2012) satisfied by RCB stars to the dereddened photometric WISE and 
2MASS data. 

Only three stars satisfying the conditions were found in the field of three globular clusters. The star in the field of Liller 1 is most 
probably a protostar. For the two other candidates, the absence of photometry in the visible range did not allow us to establish their 
nature with certainty. We further identified one dust-enshrouded star that only satisfied the first selection criteria, and used DUSTY 
to determine that it is a star of temperature 4800K enshrouded in a dusty envelope with a temperature 300 K and an opacity in the 
visible of 0.59. It is probably an Xray binary star with a dusty accretion disk. 

We found no RCB stars truly belonging to a globular cluster, thus providing a constraint on their formation mechanism. 
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1. Introduction 

R Coronae Borealis stars are a very small group of hydrogen- 
deficient and carbon-rich supergiants. About 65 are known in 
the Galaxy and 25 in the Magellanic Clouds. In our Galaxy, 
they are mostly confined to low Galactic latitudes, and are thus 
thought to be part of the bulge population. Two scenarios have 
been proposed to explain their existence, either the double de- 
generate merger of two white dwarfs, or the final helium shell 
flash of the central star in a planetary nebula (see Clayton 2012 
for a review). 

The expected number of RCB stars in the Galaxy can be as 
large as 500 (Tisserand 2012) or even 5000 (Clayton 2012). It 
thus seems feasible to increase the number of known RCB stars 
significantly, in order to better understand their evolutionary sta- 
tus. Since these stars are carbon-rich, they should be more eas- 
ily detectable from their infrared properties. Tisserand (2012) 
recently proposed eight criteria for identifying RCB stars from 
their mid-infrared WISE colors. 

Our own interest lies with the stellar populations in globu- 
lar clusters (Sharina & Davoust 2009, Sharina et al. 2010), and 
we recently discovered a CH star in the globular cluster NGC 
6426 (Sharina et al. 2012). Since the presently more likely for- 
mation scenario for RCB stars seems to be a merger of white 
dwarfs, the dense environment of globular clusters should be 
ideal for their formation. We thus decided to take advantage of 
the recently available WISE all-sky source catalogue to search 
for RCB stars in globular clusters, using the selection criteria of 
Tisserand (2012). An advantage of searching for such stars in 
globular clusters is that the distance and reddening are known a 
priori. 



2. Selection of the candidate RCB stars 

We selected the target globular clusters from the catalogue of 
Harris (1996), 2010 edition. Since RCB stars are detectable in 
the WISE bands out to 50 kpc (Tisserand 2012), we removed the 
more distant clusters. A total of 150 clusters were thus retained. 
We then extracted the WISE and 2MASS photometry from the 
WISE all-sky source catalogue, selecting all the stars within the 
tidal radius of each cluster, which was computed from the cata- 
logue of Harris (1996), 2010 edition. This produced a catalogue 
of 635989 stars. 

The magnitudes were dereddened using the standard proce- 
dure (Cardelli et al. 1989). For the four WISE bands Wl, W2, 
W3 and W4, the total extinctions R A are 0.158, 0.093, 0.087 
and 0.056, respectively (Bilir et al. 2011). The distance and 
foreground reddening for each globular cluster were taken from 
Harris (1996, 2010 edition). 

The first step was to apply the first four selection criterion 
of Tisserand (2012) to the stars. This left 2581 stars in 84 clus- 
ters. We then extracted K-band images of all the stars from the 
Interactive 2MASS Image Service and inspected them for pos- 
sible problems (artefact, nearby star, etc.) that could affect the 
photometry of each star. We considered this lengthier procedure 
safer than relying on the WISE and 2MASS flags. A total of 776 
stars and two clusters were thus removed. 

After applying the last four selection criteria of Tisserand 
(2012), we were left with only three candidate RCB stars in three 
globular clusters. Their coordinates and photometry are given in 
TableQ] None of the three stars is in the list of RCB candidates of 
Tisserand (2012, his Table 5). It is rather disturbing that all three 
candidates are situated at a very low galactic latitude (b < 0.16 



1 



Scarlett-Rose Boiardi, Sebastien Roger and Emmanuel Davoust: A search for RCB stars in globular clusters 



6 



8 



9 









□ 












w1 i 1 

w2 --X-- 
w3 

W4 : B : 


CP 


b 








□ 


B 


□ 




□ 


X- 


* 


* 


* 


* 








* 




X 


X 


X 


X 


X 


X 


X 


X 


X 


X 




















x - 


+ 


+ 


♦ 


* 


♦ 


+ 


* 


+ 


+ 


+ 


55270.4 




55270.6 




55270 


8 




55271 




55271 .2 



Fig. 1. WISE magnitudes of the candidate RCB star in Liller 1 
vs Julian Day. The dusty shell became slightly brighter in the 
redder two bands halfway through the observations. 
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Fig. 2. Spectral energy distribution of the three candidate RCB 
stars. From top to bottom : in Liller 1, 2MS-GC01 and 
GLIMPSE01. The fluxes in the four IRAC bands are also plot- 
ted. The flux units are arbitrary. 

degrees); this suggests that they are all field stars that happen 
to be in the field of a globular cluster and thus that the adopted 
distances and reddenings might be incorrect. 

3. Analysis of the three candidates 

The candidate RCB star in the field of Liller 1 has ten obser- 
vations in the WISE single-exposure catalogue. We plotted the 
four WISE magnitudes of the star against Julian Day on Fig. Q] 
which shows that the dusty shell became slightly brighter in the 
redder two bands halfway through the observations, with a hint 
of oscillatory behavior. 

An infrared color-magnitude diagram of Liller 1 has been 
obtained by Valenti et al. (2010). If we compare the magnitudes 
of the star given by Valenti et al. (2010), which date from 30 



July 2007, and those given by 2MASS, which were obtained 
on 13 August 1998, the star did not vary in H and K and be- 
came 0.5mag brighter in J between 1 198 and 2007. The K-mag 
given by DENISj, which dates from an earlier epoch than 2005, 
is comparable to the two other values. 

We obtained images in V and I of Liller 1 from the ESO 
archives, reduced them in a standard way and applied sextrac- 
tor (Bertin & Arnouts 1996) to the images. The zero-point of 
the I-band images was calibrated with photometry of 14 stars 
from DENIS. Those of the V-band were calibrated by comparing 
our color-magnitude diagram with that of Ortolani et al. (1996). 
Since the latter did not publish any catalogue, the uncertainty on 
the zero-point in V is rather large, and estimated at 0.2 mag. The 
derived magnitude in I of our star is / = 21.82, and it was not 
visible at all on the V-band image. Since one expects a V-band 
apparent magnitude in the range 17.6 - 19.6, which is easily 
reached in the V-image, this is not an RCB star if it is at the dis- 
tance of Liller 1 . In fact, the spectral energy distribution, shown 
on Fig. |2] suggests that this is a protostar (see Fischer et al. 2012 
for such a spectrum). 

We were not able to find images or photometry in the visible 
range of the two other stars. Their spectral energy distribution, 
also shown on Fig. [2] does not suggest that they could be RCB 
stars. 

4. Search for dust-enshrouded stars 

We returned to the sample of 1815 stars satisfying the first four 
criteria of Tisserand (2012) to search for dust-enshrouded stars. 
We inspected the spectral energy distribution of all the stars, and 
found only 21 which had the shape expected from a hot star 
(r* > 4000) with a warm dusty envelope (800 > T dust > 200 
K). Out of these, only one, J174104.68-534245.4 in the field of 
NGC 6397 was not within a few degrees of the Galactic plane. 

Photometry in the visible and infrared ranges and proper mo- 
tions for the star in the field of NGC 6397 were obtained using 
Vizier. It has non-zero proper motions, and thus cannot belong 
to the cluster. It has also been detected by Chandra, and is thus 
probably an X-ray binary star with a dusty torus. 

We fit a model composed of a radiation point-source and a 
dusty shell to the spectral energy distribution of the star, using 
the publicly available software DUSTY (Nenkova et al. 2000). 
We assumed a central source of temperature 4800 K and a shell 
of amorphous carbon grains with a standard MRN distribution 
of sizes (Mathis et al. 1977), visual optical depth 0.58 and tem- 
perature 300 K. The fit is shown on Fig. [3] Since we do not know 
the exact distance to the star, we plotted the fluxes corrected and 
uncorrected for extinction. 



5. Conclusion 

The negative result of our search for RCB stars in Galactic glob- 
ular clusters indicates that such stars are very rare, if not com- 
pletely absent, in these old stellar systems. This may be because 
the RCB phenomenon has a very short lifetime, is a rare event 
in the evolution of stars, or that it occurs in stars younger than 
five to seven Gyr, which is the age range of the youngest clusters 
explored in this project. 

Acknowledgements. This report summarizes the results of a summer intern- 
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Strasbourg, France. It made use of data products from the Wide-field Infrared 
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Table 1. Magnitudes and distance to the cluster center (in units of the tidal radius) for the three candidate RCB stars 



Name 


I J 


H 


K 


Wl 


W2 


W3 


W4 


cluster 


r/r, 


J173319.71-332320.1 
J180844.18-194928.4 
J184931.40-013314.7 


21.82 13.945 
17.979 
16.060 


11.753 
15.788 
13.909 


10.167 
12.876 
12.232 


9.258 
10.770 
10.972 


8.653 
10.068 
10.560 


7.832 
7.941 
8.378 


7.858 
7.211 
7.771 


Liller 1 
2MS GC01 
Glimpse 01 


0.08 
0.70 
0.79 
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Fig. 3. Spectral energy distribution of the dust-enshrouded star in 
the field of NGC 6397. The solid squares are the fluxes corrected 
for an extinction of 0.18, while the open squares are the uncor- 
rected fluxes with error bars. The solid curves show the model 
spectral energy distribution expected for a dust-enshrouded star. 
The monochromatic fluxes are normalized to the observed flux 
in the B band. 
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